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We implement an ensemble quantum counting algorithm on three NMR spectrometers with 1 H 
resonance frequencies of 500, 600 and 750 MHz. At higher frequencies, the results deviate markedly 
from naive theoretical predictions. These systematic errors can be attributed almost entirely to 
off-resonance effects, which can be substantially corrected for using fully-compensating composite 
rotation pulse sequences originally developed by Tycko. We also derive an analytic expression for 
generating such sequences with arbitrary rotation angles. 



I. INTRODUCTION 



> . 

Quantum computers are information processing devices which operate by — and exploit — the laws of quantum 
mechanics, giving them the potential to tackle problems which are intractable using classical computers [JIJ. Although 
there has been a great deal of interest in the theory of quantum computation, actually building a general purpose 
, quantum computer has proved extremely difficult. The multiparticle coherent states on which quantum computers 
rely are extremely vulnerable to the effects of errors, and so considerable effort has been expended on tackling random 
errors introduced by decoherence, including methods of error correction and fault tolerant computation [|J. 
Comparatively little effort, however, has been directed at the issue of systematic errors, arising from reproducible 
imperfections in the apparatus used to implement quantum computations. 
■ It makes sense to address systematic errors, however, because many of them can be eliminated relatively easily. It 
is convenient to visualise one spin quantum states as three-dimensional vectors, where the x, y, and z components 
are given by decomposition of the state into a linear combination of Pauli spin matrices. In this picture, unitary 
operations are interpreted as rotations of the Bloch vector on a sphere. The sensitivity of the final state to rotational 
imperfections can be much reduced by replacing single rotations with composed rotations. Ordinarily replacing one 
erroneous operation by several erroneous operations will increase the overall error, but this isn't necessarily the case 
as rotations are non-linear, and the non-linearity opens the possibility of errors cancelling one another. 

It is worth stressing the distinction between improved experimental technique, quantum error correcting codes [Q, 
and the use of composite pulses. Improved experimental technique minimises conditions which lead to data errors. 
Error correcting codes diagnose errors (however obliquely) and use this knowledge of the errors to explicitly reverse 
them. Composite pulses prevent errors in the first place by reducing the impact on the system of conditions which 
q-i cause systematic data errors, without actually eliminating these conditions. Composite pulses cannot, however, correct 
random errors, as they rely on the reproducible nature of systematic errors to correct them. 
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1 II. SYSTEMATIC ERRORS IN NMR QUANTUM COMPUTERS 

a ; 

In the last few years nuclear magnetic resonance (NMR) techniques || have been used to implement small quantum 
computers ||. Several simple quantum algorithms have been implemented using NMR, such as Deutsch's algorithm 

^| and Grover's algorithm [[io|-|l2|. The comparative ease with which such experiments have been performed reflects, 
in part, the fact that more conventional NMR experiments, used throughout the molecular sciences, themselves rely 
on the generation and manipulation of multi-spin coherent states, and so techniques for handling them are highly 
developed. 

Despite this pre-existing experimental sophistication, NMR quantum computers are nonetheless subject to system- 
atic errors. Implementing complex quantum algorithms require a large network of logic gates, which within an NMR 
implementation entails even longer cascades of pulses. In this case small systematic errors in the pulses (which can 
be largely ignored in conventional NMR experiments) accumulate and become significant. 

Interactions between individual spins and between spins and the environment are mediated by RF pulses, which are 
applications of an RF field with phase <p (in the rotating frame) [|| for some duration r. In the ideal case, the pulse 
will drive the Bloch vector through an angle 9 about an axis orthogonal to the z-axis and at an angle <fi to the x-axis. 
The rotation angle, 9, depends on the rotation rate induced by the RF field, usually written u\, and the duration of 
the pulse, r. In practice the RF field is not ideal, and this leads to two important classes of systematic errors, usually 
referred to as pulse length errors and off-resonance effects. 
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Pulse length errors arise either when the length of the pulse is set incorrectly, or (more commonly) when the RF 
field strength deviates from its assumed value, so that the rotation angle achieved deviates from its theoretical value. 
This effect is most commonly observed within NMR as a result of inhomogeneity in the applied RF field; in this case 
it is impossible for all the spins within the sample to experience the same rotation angle. 

Off resonance effects arise from the use of a single RF source to excite transitions in two or more spins which have 
different resonance frequencies. This is done for a variety of practical reasons. For example, in conventional NMR 
experiments, which seek to analyse molecular systems, there will be a large number of different transition frequencies, 
whose exact values are unknown at the start of the experiment. Thus the only practical approach is to use a single RF 
source, with sufficient power that it can excite transitions across a wide range of frequencies. In quantum computation 
experiments the transition frequencies are known beforehand, but for a variety of reasons the use of a single RF source 
for each type of atomic nucleus remains the most practical approach (different sources are used for different nuclei, 
such as 1 H and 13 C). 

Composite pulses are widely used in NMR to minimize the sensitivity of the system to pulse-length and off- 
resonance errors, by replacing direct rotations with composite rotations which are less sensitive to such effects. 
However, conventional composite pulse sequences are rarely appropriate for quantum computation because they 
usually incorporate assumptions about the initial state of the spins or introduce phase errors. These conventional 
sequences were initially derived by considering the trajectories of Bloch vectors during a pulse, and this approach is 
only useful if the starting point of the Bloch vector is known. 

QC requires fully-compensating (type A) composite pulse sequences [[l3|. Fortuitously, a number of such sequences 
were developed fifteen years ago by Tycko Q. While these pulse sequences do not offer quite the same degree of 
compensation as is found with some more conventional sequences, the compensation which does occur is effective 
whatever the initial position of the Bloch vector. As such sequences are rarely (if ever) required within conventional 
NMR experiments, these sequences have previously found no experimental use. They are, however, ideally suited to 
quantum computations. 



III. EXAMPLE: NMR QUANTUM COUNTING 



Systematic errors are an issue with many types of quantum algorithms, but we will focus on a counting algorithm 
which can involve particularly long pulse trains. Counting the number of items matching a search criterion is a 
well-known computer science problem, and an efficient quantum counting algorithm based on a variation of Grover's 



quantum search has been implemented using NMR techniques 1 15 



Consider a match function f(x) which maps n-bit binary strings to a single output bit. In general there will be 
N = 2" possible input values, of which k will yield a match, /(x)=l. The object of the quantum counting algorithm 
is to estimate k by estimating an eigenvalue of the Grover iterate G = HUqH^Uj, where H represents an n-bit 

Hadamard transform, Uq maps |0) to — |0), and Uj maps \x) to (— 1)/( ;E )+ 1 \x). The algorithm involves applying G 
repeatedly and observing the variation in signal intensity. The signal intensity is modulated in r, the number of 
applications of G, with frequency proportional to k, the number of matches. A more detailed explanation of this 
algorithm is found in reference [J15[ . 
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FIG. 1. Experimental results from a 500 MHz NMR quantum computer for two of the four possible functions /: (a) foo 
(no matches); (b) fn (both match). The observed signal intensity is plotted as a function of r, the number of times the 
controlled-G operator is applied. Intensities are normalised relative to the case of r = 0. The solid lines are exponentially 
damped cosinusoids with the appropriate theoretical frequencies, and are plotted to guide the eye. 
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Figure |l] shows experimental results for a search over a one qubit search space, using a two-qubit quantum computer 
based on the two 1 H nuclei of cytosine in solution in D2O ||. The NMR spectrometer used a 1 H operating frequency 
of 500 MHz. On a one qubit search space, there are three possible results corresponding to four possible functions /: 
no inputs match (in which case / = /oo)> the first input matches (/ = /01), the second input matches (/ — /10), or 
both inputs match (/ = /n). For simplicity, only the cases /oo and fn are shown. 

The observed signal loss with r could arise from a number of sources. Clearly one possibility is the effects of 
decoherence, but the observed decay rate is too rapid to be so simply explained, and patterns in the signal loss (this 
is especially clear for /oo, where data points lie alternately above and below a smooth curve) suggest a more complex 
explanation. Repeating the experiments on spectrometers with higher operating frequencies markedly degrades the 
results, as shown in figure which shows data from computations implemented on 600 and 750 MHz spectrometers. 
Once more, only the /oo and fn cases are shown. These results clearly demonstrate that the use of expensive high 
field magnets does not always give better results. An unwanted beat frequency and high frequency chatter have been 
introduced, and the signal decays more rapidly. These effects are more severe at 750 MHz than at 600 MHz. 
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FIG. 2. Experimental results from a 600 MHz ((a) and (b)) and 750 MHz ((c) and (d)) NMR quantum computer for / 00 ((a) 
and (c)) and fu ((b) and (d) ). The observed signal intensity is plotted as a function of r, the number of times the controlled-G 
operator is applied. Intensities are normalized relative to the case of r = 0. 

The observed field dependence strongly suggests that these errors arise from off-resonance effects. These arise 
because only a single RF transmitter is used to excite both 1 H nuclei, and it cannot be on-resonance with both sets 
of transitions. Instead the transmitter is placed halfway between the two resonance frequencies, resulting in equal 
and opposite off resonance terms. Two factors are responsible for the poor results at higher frequencies. Firstly, the 
frequency offsets, ±<5i/, are proportional to the resonance frequency, and thus are larger at high frequencies. Secondly, 
the RF field strength is typically weaker at higher frequencies. These combine to make the off-resonance effects 
much more serious. This diagnosis can be confirmed by including off-resonance effects in simulations of the NMR 
experiment. After this correction is made, the theoretical results, shown in figure [|, agree strikingly well with the 
experimental results. Even the high frequency 'noise' on the peaks is reproduced. 
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FIG. 3. Theoretical prediction for 750 MHz spectra once off-resonance effects are included. An (arbitrary) exponential 
damping has been applied. 

Since our simulations confirm that the major errors in NMR quantum counting arise from off-resonance effects in 
the 90° pulses, we chose to use Tycko's 385° 320°_ y 25° composite 90° sequence, which offers good compensation for 
off-resonance errors while remaining fairly insensitive to pulse-length errors. Substituting this sequence for the 90° 
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pulses produced a significant improvement in the experimental results at 750 MHz, as shown in figure The low 
frequency beating almost disappears, and the high frequency noise is much reduced. 



0.5 - \ 



16 32 48 64 80 96 112 128 16 32 48 64 80 96 112 128 

(a) (b) 

FIG. 4. Experimental results from the 750 MHz quantum computer when 90°, pulses are replaced by 385^320-1, 25 y pulse 
sequences. 

When using composite pulses, there is a trade-off between the cancellation of errors and the extra manipulation 
of the system necessary to induce this cancellation. While there exist sequences which in theory compensate for an 
extraordinary range of errors, these sequences involve prohibitively long cascades of pulses. In practice, three-pulse 
sequences seem to provide a good balance between insensitivity and simplicity. The hazards of over-manipulation can 
be seen when the counting experiment is repeated using composite 180° pulses as well as 90° ones, shown in figure 
IL The composite 180° pulse used was 90° 225°_ a , 315° |ll|. Although simulations predict that the results should be 
slightly better when both 180° and 90° compensated pulses are used, the results for /n are in fact slightly worse. 
The errors introduced by the extra manipulation outweigh the small gain achieved by the use of more sophisticated 
pulse sequences. 
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FIG. 5. Experimental results on the 750 MHz spectrometer when 90°, pulses are replaced by 385°, 320 c L y 25°, pulse sequences 
and 180°, pulses are replaced by 90°, 225°.^ 315°, pulse sequences. 

IV. COMPOSITE PULSES FOR ARBITRARY ROTATION ANGLES 

Discussion of composite pulses in the NMR literature has focussed almost exclusively on 90° and 180° pulses, for 
the simple reason that these are the rotation angles most commonly used in conventional pulse sequences. While 
these composite pulses can be very useful for quantum computation, as shown above, it would be preferable to have 
composite pulses for other rotation angles as well. For example, a 45° pulse could be used in the implementation of 
a true Hadamard gate Jl7t , which might be applied many times in a single network. 

We use the method of coherent averaging to generate an analytic expression for composite pulses similar to Tycko's 
90° pulse but with arbitrary rotation angle. The basic idea of the method is to split the Hamiltonian of the system 
into intended and error components, 

H(t)=H (t) + V, (1) 
where the "ideal" Hamiltonian (written in the rotating frame and using product operator notation) is 
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H (t) — uii (I x cos fa+ I y sin fa (2) 

and the error term is 

V = 8ljI z (3) 

{5w is the resonance offset angular frequency). We then seek to minimise the propagator of the error component by 
expanding it as a power series Knowing V , we can rewrite the propagator as a product, 

U{t) = U (t)U v (t) (4) 

= Texp(-i/ dtH (t))Texp(-i dtV{t)) (5) 



where T is the Dyson time-ordering operator and 

V(t) = UoityWUoit). (6) 
Since Uy varies rather slowly with time, it can be expanded as a power series in V: 

U v (t) = exp f-ir (V (0) + V (1) + ••■)) (7) 



using the Magnus expansion 



V m =- f dt V(t) (8) 

T In 



V™ = -£ I dt, I dt: 



ti 



Vit^^ih) (9) 



etc. Complete expressions for the higher order terms are given in |18j,|l9|]. Our objective is to find a pulse sequence 
which satisfies the requirement that it perform an ideal rotation under ideal conditions, 

m 

TjlUiin) = U (t) (10) 

i=l 

while minimising V^°\ We do not consider higher order terms as using more equations would require more variables in 
the solution and we wish to restrict the length of our composite sequence to three pulses. Here Uq(t) is the propagator 
of an ideal uncomposed pulse and the Ui(n) are the constituent propagators of the ideal composite pulse. Ordinarily, 
a numerical search for a solution must be conducted for each set of target values, 9 and </>, but an analytic solution 
exists if the axes of rotation are taken to be <f>i — <f>, fa = 7r + fa and fa = <fi. While this is nominally a loss of 
generality, the solutions we find are perfectly compensated to first order, and so no more general solution would be 
better to first order. The solutions also have reasonable properties with respect to RF field inhomogeneity. 
With these phases, it follows trivially from equation (|l^) that 

e 2 = e 1 + e 3 -e (11) 

Working out with these values of (f> and 02 and setting Vy ^ and to zero (it happens that Vx°^ is already 
zero) gives a restriction on 9$ , 

9 a = ±9 1 + 2nir n = 0, ±1, ±2, . . . (12) 
Taking the positive result, 9% = 6\ + 2n7r, and back-substituting gives two solution families, 



f (l-cosfl) 2 ±sinflV(l-cosfl)(7 + cos^ , - m > 

i\ = iarccos ; — + 2ri7r (13 

\ 4(1-008 6*) / 



sin 9(1 + cos ( 



^(l + cos(9) (7 + cos (9)) 



iarcsin ^ — + 2mr (14) 

1 4(l + cos6>) 1 
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These solutions are equal in the region > 8 target > n for the positive root and a consistent choice of signs. (This is 
easily checked by working out the Taylor expansions about 0, and less easily checked by working out the difference 
between equations ( |l3| ) and (|TJ) using an identity for arcsina — arccos/3.) 

Although we have not explicitly considered the effects of RF field inhomogeneity when determining the 4>i and 9\ , 
it happens they are trivial for the rotation axes we have chosen. Using V r f = 5u>\ (I x cos <f> (t) + I y sm. <f> (t)) , where 
Suji is a measure of the field inhomogeneity, for V and working out with the <j>i specified above, we get 

V^=6u x (15) 
= (16) 

uj 0) =o (17) 



Therefore, equations (13) and ( |14|) can be used exclusively to determine 9\. We follow Tycko and choose n — 1 so 



that 2tt < 0i < 4tt. The final rotation angles are 



/ (1 - cos 6Y + J 1 - cos 6>) 7 + cos 9) sin 9 \ 

8i = arccos + 2tt 

I 4(l-cos0) J 

9-2 = 0\ + 03 — 9 

03 = 01 - 2tt (18) 
Within this family, the quality of the pulses is remarkably insensitive to to actual details of the pulse sequence. As 



long as equation (|ll|) is strictly satisfied, any sequence with angles vaguely like those given by equation (18) will 
produce broad-band rotations. 

Figure ^| shows fidelity plots for Tycko's original 90° composite pulse, and a new pulse derived using our analytic 
method. For comparison purposes, fidelity plots for the equivalent uncompensated pulses are included as well. The 
fidelity was evaluated using quaternions as described in reference Q . 

The operation of the sequence whose fidelity is plotted in Figure ||(d) is illustrated in Figure fj]. The resonance 
offset was chosen to be fairly large so that the mechanism of the composite pulse is clear. 

The composite pulses we have derived have good properties in both quantum computation and spectroscopic 
contexts, but it may be possible to produce pulses with properties that are even more favourable for quantum 
computation if we stop using spectroscopy as our implicit model. Consider equation (|l^) again. Although it seems 
like an obvious condition that the behaviour of the pulse should be ideal in the absence of errors, this is misleading. 
We would never want to use a composite pulse in a situation where all pulses were on-resonance, and so there is no 
a priori reason to design pulses which perform well in the event of this non-existent ideality. The resonance offsets 
for quantum computation will always be known exactly, so we would like to be able to design pulse sequences which 
perform best as these known offsets rather than at the zero offset. This is a problem we are currently considering. 



V. CONCLUSIONS 



Composite pulses show great promise for reducing data errors in NMR quantum computers. More generally, any 
implementation of a quantum computer must be concerned on some level with rotations on the Bloch sphere, and 
so composite pulse techniques may find very broad application in quantum computing. Composite pulses are not, 
however, a panacea, and some caution must be exercised in their use. 
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FIG. 6. Numerical evaluation of some composite pulses. The :r-axis represents radio-frequency inhomogeneity and the 
y-axis off-resonance effects. The sequences are (a) An uncompensated 90J pulse (b) Tycko's 90J pulse, 385^ 320^ 25J (c) 
an uncompensated 60J pulse (d) a new 60^ composite pulse, 375^ 331^ 15^ In all of these sequence the axes of rotation are 
4>i = 0, <j>2 — 7r, cj>3 — 0. Solid contour lines are spaced at 15% intervals, and dashed contour lines are spaced at 1% intervals, 
beginning at 95%. 




(a) (b) 

FIG. 7. Magnetization trajectories on the Bloch sphere for a 60° uncompensated pulse and Tycko-type composite pulse for a 
system with a fairly large resonance offset. The black dot shows the final position on the Bloch sphere, and the pale ring shows 
the target position. The three rotations within the Tycko-type sequences are shown in different colours, (a) An uncompensated 
rotation of a spin system with a resonance offset of 2 ppm in a 750 MHz magnet with an RF field strength of 10 KHz. (b) A 
composite rotation in the same system. 



7 



[1] P. W. Shor. SI AM Review, 41:303, 1999. 
[2] P. W. Shor. Physical Review A, 52:2493, 1995. 
[3] A. Steane. Physical Review Letters, 77:793, 1996. 
[4] A. Steane. Nature, 399:124, 1999. 

[5] R. R. Ernst, G. Bodenhausen, and A. Wokaun. Principles of Nuclear Magnetic Resonance in One and Two Dimensions. 
Oxford, 1997. 

[6] D. G. Cory, A. F. Fahmy, and T. F. Havel. Nuclear magnetic resonance spectroscopy: an experimentally accessible paradigm 
for quantum computing. In T. Toffoli, M. Biafore, and J. Leao, editors, PhysComp96, page 87. New England Complex 
Systems Institute, 1996. 

[7] D. Deutsch. In R. Penrose and C. J. Isham, editors, Quantum Concepts in Space and Time. Clarendon, Oxford, 1986. 
[8] J. A. Jones and M. Mosca. J. Chem. Phys., 109:1648, 1998. 

[9] I. L. Chuang, L. M. K. Vandersypen, X. Zhou, D. W. Leung, , and S. Lloyd. Nature, 393:143, 1998. 
[10] L. K. Grover. Physical Review Letters, 79:325, 1997. 

[11] I. L. Chuang, N. Gershenfeld, and M. Kubinec. Physical Review Letters, 80:3408, 1998. 

[12] J. A. Jones, M. Mosca, and R. H. Hansen. Nature, 393:344, 1998. 

[13] M. H. Levitt. Progress in Nuclear Magnetic Resonance Spectroscopy, 18:61, 1986. 

[14] R. Tycko. Physical Review Letters, 51:775, 1983. 

[15] J. A. Jones and M. Mosca. Physical Review Letters, 83:1050, 1999. 

[16] Z. Starcuk and V. Sklenar. Journal of Magnetic Resonance, 62:113, 1985. 

[17] J. A. Jones, R. H. Hansen and M. Mosca, Journal of Magnetic Resonance, 135:353, 1998. 

[18] U. Haeberlen and J. S. Waugh. Physical Review, 175:453, 1968. 

[19] U. Haeberlen. High Resolution NMR in Solids. Academic Press, 1968. 



8 



